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Abstract 



The technique of reduction of order developed by Nucci (J Math 
Phys 37 (1996) 1772-1775) is used to produce nonlocal symmetries 
additional to those reported by Krause (J Math Phys 35 (1994) 5734- 
5748) in his study of the complete symmetry group of the Kepler Prob- 
lem. The technique is shown to be applicable to related problems 
containing a drag term which have been used to model the motion 
of low altitude satellites in the Earth's atmosphere and further gener- 
alisations. A consequence of the application of this technique is the 
demonstration of the group theoretical relationship between the simple 
' harmonic oscillator and the Kepler and related problems. 

I-H . 

1 Introduction 

'a ■ 

^1 In a paper of a few years ago Krause |15] introduced a new concept into 

the study of the symmetries of ordinary differential equations. He called 
this a complete symmetry group and defined it by adding two properties 
5^ I to the definition of a Lie symmetry group. These were that the manifold 

of solutions is an homogeneous space of the group and the group is specific 
to the system, ie no other system admits it. This definition required the 
introduction of a new type of symmetry defined by 



1 ^ 

Y= ^{t,xi,...,XN)dt dt + '^'ni{t,xi,...,XN)dx^. (1.1) 

1=1 
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This definition of a symmetry differs from that of a Lie point symmetry due 
to the presence of the integral as the coefficient function of dt- 

As an illustration of the concept of a complete symmetry group Krause 



used the Kepler problem and obtained three symmetries of the type of (1.1)- 
He claimed that these three symmetries could not be obtained by means of 
the standard Lie point symmetry analysis. Naturally it was not long be- 
fore this claim was shown by Nucci ||25| to be incorrect in the case of an 
autonomous system. In the case of the Kepler problem, an autonomous sys- 
tem, one of the dependent variables can be taken to be the new independent 
variable and the order of the system be reduced by one. An analysis of the 
reduced system for Lie point symmetries leads to results different from the 
analysis of the original system. In particular the three additional nonlo- 
cal symmetries obtained by Krause followed from point symmetries of the 
reduced system. 

One of the fundamental problems of mechanics is that of the Kepler 
problem which describes the interaction of two point particles with an inverse 
square law of attraction. It is well-known that this problem possesses the 
first integrals of the conservation of the scalar energy, the vector of angular 
momentum and vectors in the plane of the orbit known as Hamilton's vector 
[10 1 and the Laplace- Runge-Lenz vector |ll], ||, |l^, 27, 22]. The invariance 



Lie algebra of the first integrals under the operation of taking the Poisson 
Bracket is so(4) (in the case of negative energy) and the Lie algebra of the 
five Lie point symmetries of the equation of motion A2 so(3). The algebra 
of the complete symmetry group has not been given. The elements of the 
five-dimensional algebra are 

Xl = dt X3 = X2dx3 - X3dx2 

X2 = tdt + Irdr Xi = xsdxi - xid^^ (1.2) 
X^ = xid^^ - X2dx^. 

The additional three nonlocal symmetries provided by Krause are 
¥1 = 2(^1 xidt^ dt + xirdr 
¥2 = 2(^1 X2d?j dt + X2rdr 

Yi = 2 (^1 xsd^j dt + xsrdr (1.3) 
in which = Xi + x"^ + Xn. 
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There have been several other systems, generahsations of the Kepler prob- 
lem, which have been shown to have a similar set of conserved quantities 
[^, ^, ^ . Just as the Laplace- Runge-Lenz vector provides a direct route 
to the equation of the orbit of the classical Kepler problem, the correspond- 
ing vectors of the generalised Kepler problems provide the same direct route 
to the equations of their orbits. The Lie point symmetry associated with the 
Laplace-Runge-Lenz vector of the classical Kepler problem is the rescaling 
symmetry, X2. The generalisations of the Laplace-Runge-Lenz vector do 
not always have such an associated Lie point symmetry. In the case of the 
equation of motion 

- f|- + V + = 0, (L4) 



2g 2r 

which is a variation of the model proposed by Danby |2^, |l9| for the 

motion of a satellite in a low altitude orbit subject to atmospheric drag, viz 

ar ur , ^ 

r + 7^ + 73=0' (1-5) 

and for which the generalisation of the Laplace-Runge-Lenz vector is 

J = — + -^, (1.0) 

where L is the magnitude of the angular momentum, L, and ^ is a constant 
of the motion defined through 

1 

L = A[gr^y, (1.7) 

Pillay et al [^] showed that, instead of the Lie point symmetry for the 
classical Kepler problem, X2, the nonlocal symmetry 



^dt 

g 



dt + rdr {II 



was the corresponding associated Lie symmetry. 

In this paper we intend to demonstrate the existence of far more nonlocal 
symmetries for the classical Kepler problem than were reported by Krause. 
We derive them as Lie point symmetries of a reduced system using the 
method of Nucci ||2^ . By deriving these additional nonlocal symmetries in 
this way we are able to say something definite about the expanded symmetry 
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group. The essence of the method of Nucci is to reduce the order of the 
system by usmg the symmetry, Xi, which is a statement of the autonomy 
of the system. We recall that in a reduction of order the symmetry Z which 
does not have the property [Xi, Z] = \Xi becomes an exponential nonlocal 
symmetry By an exponential nonlocal symmetry we mean one of the 
form 



G = exp 



fdt 



{rdt + riA;) (1.9) 



in which without a knowledge of the solution of the differential equation / is 
not an exact derivative and r and the r]i are functions only of t, the Xi and 
their derivatives. In this case we have nonlocal symmetries becoming local 
on the reduction of order. Consequently we know that the symmetries of 
the reduced equation quite possibly have zero Lie Bracket with Xi and this 
will enable us to construct the algebra. Further we shall show that these 
considerations which are applicable to the classical Kepler problem can be 



extended to the generalisations such as the one in (1.4). 

In the next section we review the reduction procedure of Nucci |25| and 
in Section 3 we apply it to the classical Kepler problem and see that there 
is a certain delicacy in the choice of the new independent variable. For the 
sake of simplicity we work in two dimensions. In Section 4 we make some 
observations about these symmetries, the route to further simplification and 
the algebra. In Section 5 we obtain the results for Danby problem Q, in 
Section 6 those for the generalised problem represented by (|1.4| ), in Section 
7 the symmetries for another generalisation in which the force is not only 
not central but is also angle dependent and in Section 8 we present our 
conclusions and make some pertinent observations about going to the full 
three dimensions. 

2 The method of reduction of order 

Consider the system of second order ordinary differential equations given 
by 

= /,(x,i;), i = l,N (2.1) 

in which t is the independent variable and = 1,A^ the dependent 
variables. These equations may be considered as equations from Newtonian 
mechanics, which was Krause's approach, but there is no necessity for that 
to be the case. There is also no necessity for the dependent variables to rep- 
resent cartesian coordinates. Indeed there is no need for the system to be of 
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the second order. It just so happens that many of the equations which arise 
in practice have their origins in Newton's Second Law and so are second 
order equations. There is no requirement that the system be autonomous. 
In the case of a nonautonomous system we can apply the standard proce- 
dure of introducing a new variable xtv+i = t and an additional first order 
equation xat+i = 1 so that the system becomes formally autonomous. In 
our discussion we confine our attention to autonomous systems. We reduce 



the system (2.1) to a 2A^-dimensional first order system by means of the 



change of variables 



Wl = Xl Wn+1 = Xl 

W2 = X2 Wn+2 = X2 

WN-1 = XN-1 W2N-1 = XN-1 

Wn = X]\[ W2N = Xn 



(2.2) 



so that the system ( |2.lD becomes 

Wi = gi{x,w), i = l,2iV, (2.3) 

where gi = wm+i for i = \,N and gi = fi for i = N + I, 2N. 

In the first step of the reduction of the original system (2.1) we simply 



follow the conventional method used to reduce a higher order system to a 
first order system. Any optimisation is performed in the further selection of 
the final variables. This selection may be motivated by the existence of a 
known first integral, such as angular momentum, or some specific symmetry 



in the original system (2.1) 



We choose one of the variables Wi to be the new independent variable 
y. For the purpose of the development here we can make the identification 
Wn = y- By taking the quotients of the first order equations of the remaining 
members of the set ( |2.2D with ( |2.3| N) we obtain the (2A^ — l)-dimensional 
system 

^""^ -• = l,...,A^-l,iV + l,...,2iV. (2.4) 



dwAT gN W2N 

We do not attempt to calculate the Lie point symmetries of the system 
( |2.4D because the Lie point symmetries of a first order system are gener- 
alised symmetries and one has to impose some Ansatz on the form of the 
symmetry. Rather we select n < N — 1 the variables to be the new depen- 
dent variables and rewrite the system (|2.4D as a system of n second order 
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equations plus 2{N — n) — 1 first order equations. The selection of the new 
dependent variables is dictated by a number of considerations. The first and 
foremost is that we must be able to eliminate the unwanted variables from 
the system (2^). After this condition has been satisfied we may look to 
seek variables which reflect some symmetry of the system, for example an 
ignorable coordinate such as the azimuthal angle in a central force problem. 

After the symmetries have been calculated, they can now be translated 
back to symmetries of the original system as follows. Suppose that the 
symmetry in the original variables is given by 

G = Tdt + md^,. (2.5) 

The symmetry G is first extended and then rewritten in terms of the new 
coordinates as follows 

= Tdt + r/i^.^ + {f]i - Xif)d±^ 
= Tdt + Cidw, 

= ady + ^idu„ (2.6) 

where in the first line the summation is from 1 to A^, in the second from 1 to 
N — 1 and A + 1 to 2A^ and in the third over the number of dependent vari- 
ables Ui (the number cannot be fixed in advance without a knowledge of the 
specific system); Q = rji for i = 1, N — 1 and Ci = 'ni~ 9i'^ for i = A + 1, 2A; 
c" = ?77v; = Cjdui/dwj. The only way that r appears in the symmetries of 
the reduced system is through its derivative with respect to time. If the non- 
locality in the original system occurs as a simple integral in r of a function 
of the original dependent variables, Xj, this will be passed to the reduced 
system as a function of the new variables. When the point symmetries of 
the reduced system are computed, the form which the symmetries take in 
the original system can be determined from ( |2.6| ). Since r is determined 
as its derivative with respect to time, the symmetry of the original system 
must necessarily be nonlocal unless the derivative is an exact differential. 
We note in passing that there is no inherent restriction on the nature of the 
symmetries. They could equally be contact or generalised symmetries and 
the same considerations would apply. The only requirement is that r be a 
simple integral, not that the integrand be a point function. However, for 
the purposes of this paper we confine our attention to point symmetries and 
the integrand in r to a point function. 

In addition to the nonlocal symmetries which may be collected by this 
procedure the reduced system will have as point symmetries those symme- 
tries of the original system which have the correct Lie bracket with the 
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symmetry dt which is at the basis of the reduction of order outhned above. 
Thus, if G is a symmetry of the original system and 



[G,dt]=Xdu (2.7) 

G wiU, when expressed in the appropriate coordinates, be a point symmetry 
of the reduced system whereas, if 

[G,dt]^Xdt, (2.8) 

G wiU not be a point symmetry of the original system but an exponential 
nonlocal symmetry [Q]. Consequently there is the potential for a loss of 
symmetry in the reduction process just as there is the hope of an increase 
in the total number of symmetries, both point and nonlocal, known for the 
original system. 

3 Lie point symmetries of the reduced Kepler prob- 
lem 

The Lagrangian for the two-dimensional Kepler problem is 

L = ^(^r^ + rW)+^ (3.1) 
in plane polar coordinates and the two equations of motion are 

r-r^2 = -ii (3.2) 
re + 2re = 0. (3.3) 

We introduce the new variables and their time derivatives 
wi = r wi = 

W2 = 7112 = Wi 

W3 = r ws = wiwl - (3.4) 



w 



■ . 2W3W4: 

W4, = t) Wi = 



Wl 

In accordance with the development in the previous section we select W2 



to be the new independent variable y. The left side of (3.4) leads to the 
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reduced system 



dwi 


_ W3 


dy 




dW3 




dy 


= WIW4 




dW4 




dy 


Wl 



(3.5) 
(3.6) 

W1W4 

(3.7) 

^Prom ( p.5D we have w^ = w^w'i, where the prime denotes differentiation 
with respect to the new independent variable, y, and we replace (13.71) by 

dW4 2w4w'i 

— — . (3.»j 

ay wi 

In (^) we replace W3 by W4w[ to obtain 

„ W4W'1 fl 

W4W1 = W1W4 K — . (3-9) 

Wl wfw4 

By our replacement of ^3 we have not precisely decided that the variables 
ui and U2 are to be wi and W4. If we do make this identification, we obtain 
a system of two equations, one of the second order and one of the first order, 
viz 

12 

u'( = 2^+ui--^ (3.10) 

Ul ufU2 

Un = —2, . ("J-llj 

Ul 



We observe that (3.11) is trivially integrated to give u\u2 is a constant. (This 
is a consequence, naturally, of the symmetry dg of the original system (|3.2| ) 
and (|3.3D which is a reflection of the fact that 9 is an ignorable coordinate.) 



Consequently we may just as well define our new variables to be 

Ul = Wl and U2 = u\u2 (3.12) 
so that the system of equations we are to consider is 

u'i = 2^+ui-^ 

Ul U2 

u'2 = 0. (3.13) 
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Hereafter we drop the tilde. 

We calculate the Lie point symmetries of the system ( [3.13 ) using the 



well-known interactive program developed by Nucci and obtain the sym- 
metries 

Xl = dy 

X2 = 2uidu^ + U2(?„2 

X3 = Ul (^fiui - ufj du-, 

X4 = u\ cos ydu^ 

Xq = {^lUi — M2) COS ydy + Ul (2^ui — u^^ s'mydui 

X-j = — [jiui — sinydy + ui — u^^ cosydm 

Xg = U2 cos 2ydy — Ul (^^ui — u^^ sin 2ydui 

Xg = U2 sin2ydy + ui (^fiui — U2^ cos2ydui- (3-14) 

Some of these symmetries are readily identified, but not many of them. 
Clearly Xi represents the rotational invariance of the system and consti- 
tutes the subalgebra, so{2), of the original system of equations. In X2 we 
recognise the rescaling symmetry closely associated with the Laplace-Runge- 
Lenz vector. 

The other symmetries are not so easy to identify without making some 
calculation. To obtain the form of the symmetries in the original coordi- 



nates we must make use of (2.6). If we write the symmetry in the original 
coordinates as 

G = Tdt + 7jdr + Cde, (3.15) 
in terms of the new variables, the symmetry has the form 



G = Cdy+ rjdui + 



c + 



2?? 

r 



r^du,, (3.16) 



where the coefficient functions are all expressed in terms of the original 
variables. The calculations are not particularly interesting and we simply 
list the symmetries. They are 

Xi = de 

X2 = 3tdt + 2rdr 
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^3 
Xi 
^5 
Xq 

Xj 

Xs 
Xg 



fi J r'^dt - L^t 
2 y r cos Odtdt + cos Odr 
2 / r sin Odtdt + sin Odr 



dt + fir [r"^ — L"^) (9, 



3^r — L^j sin 9 + ^- cos i 
+ cos 6* (^/ir — L^^ 



dtOt + r sin (2/ir - L^) 

3^r — L^^ cos — fi^ sin dtSt + r cos ^2;ur — L^^ 5^ 
— sin (^fj,r — L^^ 

2/i y r sin 2edtdt + r (^fir - L^) sin 2edr - cos 26'(9e 
2/i y" r cos 26'dt(9t + r (//r - L^) cos 26*5^ + sin 2ede (3.17) 

so that we see that we lost only the time translation symmetry in the re- 
duction of order and so gained six additional symmetries. 

4 The Lie point symmetries of the reduced Kepler 
problem: further considerations 

Further interpretation of the Lie point symmetries of the reduced generalised 
Kepler problem is facilitated by the redefinition of the symmetries given in 
( |3.14 ). We do not alter the definitions of the first three symmetries, but we 
shall include them in this new listing for the sake of completeness. We now 
have the set of symmetries 



^1 

X2 
X3 

Xi± 

x^± 

Xq± 



dy 

2uidui + U2du2 

Xi ± iX5 = e^'yuldu. 



Xe ^ iXr = e^'y 
Xs ± iXg = e^^'y 



{jxui - ufj dy =F ui (2fxui - U2 ) <9„i 



(4.1) 
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We may search for the first integrals/invariants associated with each of these 
symmetries in the usual way. By way of concrete example we take X4_|_ . The 
invariants of the first extension of X4_(_, viz 



X 



[1] 



are found from the associated Lagrange's system 



dy 




dui 



dn'i 



dU2 _ 
2uiu[ + iui 



and are 



a = y, P = U2 and 7 



Ul + iui 



(4.2) 



(4.3) 



(4.4) 



where the first two are by inspection and the third comes from the solution 
of the second and fourth of ( [4.3| ). The integral/invariant is a function of 
these three arguments and is found by demanding that the total derivative 
with respect to y be zero when the differential equations ( 3.13 ) are taken 
into account. We obtain the associated Lagrange's system 



da 

T 



dp 

T 



d7 



/32 



(4.5) 



which gives /? as one of the characteristics and 



Ul + iui ifj, 



ut, 



(4.6) 



as the second characteristic. Both /? and to are first integrals/invariants 
of the system ( 3.13 ). Naturally we recognise the former as the angular 
momentum. 

If we perform the same calculation with X^^, we obtain a similar result, 



VIZ 



OJ- 



-ly 



lUl 



7/2 

Ul 



III 

+ 4 



(4.7) 



The two can be combined into one convenient expression given by 

J± = e^^y {v'^iv), v = ^l- (4.8) 



where we have made use of the constancy of U2 to write J± = U2i^±. For the 
system ( p. 13 ) J± are two invariants and for the original system, ( |3.1| ) and 
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( |3.2D , the two components of the first integral known as the Laplace- Runge- 
Lenz vector, written in complex form. 



In (4.8) we introduced a new variable v. If we use this variable instead 
of ui in the system ( [j.l3 ), we obtain the system 



v" + v = 
= 



u'2 



(4.9) 



so that the natural variable which arises from the invariants of X4± is a 
variable which further simplifies the reduced system. This is an interesting 
phenomenon for we are obtaining natural variables as we progress through 
the process of determining the symmetries of the reduced equation. It be- 
hooves us to rewrite the symmetries in terms of this new variable. We find 
that 



X2 
X3 
Xi± 
Xq± 
Xs± 



du2 

dy 

vdy mod(ii2) 
e^'^d^ mod(n^) 
e^'^'y [dy ± ivd^] modiul 
vdy ± iv^d^ mod{u 



''2J- 



(4.10) 



which is certainly a simpler appearance. 

In the simpler form presented in (4.11) we see that the calculation of the 
integrals/invariants is simpler. For example, if we take Xq±, the first set of 
characteristics comes from the solutions of the associated Lagrange's system 



dy 
1 



dv 
iv 



dv' 



-IV 



2v 



(4.11) 



The characteristics are 



a 



ve 



and p 



' -■- • 2 

VV ^ IV . 



(4.12) 



The condition for the function to be an invariant (in this case not a first 
integral of the reduced system since the independent variable is explicitly 
present) is that these two characteristics satisfy the first order equation 



da 
a 



dp 



(4.13) 



12 



whence the invariants are given by 



= ^ = (v' Tw)e^'y (4.14) 
a 

which are, of course, the two components of the Laplace-Runge-Lenz vector. 
We note that there is also a first integral associated with each of Xq± and 
this is the angular momentum which, in these coordinates, is an ignorable 
coordinate. 

For the sake of completeness we list the first integrals/invariants associ- 
ated with the symmetries listed in ( 4.101 ). They are given in the same order 



and with the same subscripts as the symmetries are listed. 

Gi Ii = L l2 = J+J- = 2L^E + fi^ 

G2 Ii = L I2 = -J- 

= ^ ^' = 1^ (4.15) 

G4± Ii = L l2± = J± 

Csi Ii = L l2± = — 
Gq± Ii = L l2± = J±. 

We note that in the second integral associated with Gi we have expanded the 
product into the standard expression relating the square of the magnitude of 
the Laplace-Runge-Lenz vector with the energy and the angular momentum. 
For the second integral of G5-1- we could have equally written I2 = J^/ J^, 
but we chose to maintain the pattern of it. 



5 The Kepler problem with drag 

In the introduction we referred to the model proposed by Danby for the 
motion of a low altitude satellite subjected to a resistive force due to the 
Earth's atmosphere described by the equation of motion 

r+^ + ^ = 0, (5.1) 

where a and ^ are constants. Since the direction of the angular momentum 
is a constant, we may analyse the problem in two dimensions using plane 



13 



polar coordinates, {r,6). The two equations of motion are 



ar 



+ ^ + ^ = 



rO + 2r(9 + 



0. 



(5.2) 
(5.3) 



We introduce the new variables and their time derivatives 



wi= r 


Wl 


= W3 


W2 = 9 


W2 


= 


11)3 = 7- 


W3 


2 ' 


W4, = 9 


W4 





aw3 



(5.4) 



Wl 



We again select W2 to be the new independent variable y. The right side of 
( |5.4D becomes 





_ W3 




dy 






dW3 




aw3 


dy 


= W1W4 — 






W1W4 


dW4 


2W3 


aw4 


dy 


Wl 


wf 



wfwi 



(5.5) 
(5.6) 
(5.7) 



In this case the choice of 
the two-dimensional system 



W1W4 + w'lw'i 



to eliminate ^3 is obvious and we obtain 



W1W4 



At 



wfwi 



2w'iW4 a 



2 ■ 



Wl Wl 

In the case of ( |5.9[ ) we can easily manipulate it to obtain 

-ay + P, 



{wlw4^ 



a 



W1W4 



(5.8) 
(5.9) 

(5.10) 



where /3 is a constant of integration, which indicates that the angular mo- 
mentum is not conserved. We have a choice of defining the new variable 
U2 as either wfwi, which is not conserved but is a convenient variable for 
manipulations, or wlw^ + ay, which is conserved but is not a convenient 
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variable for manipulation. For the present we make the former choice. In 
equation (^) we make use of ( 5.10 ) to eliminate W4 and 1^4 to obtain 



1 

Wl 



1 

Wl 



iay+(3) + 



1 

Wl 



We introduce the second new variable as 

ay + P 



Ul 



Wl 



{ay+[3f+li = Q. (5.11) 



(5.12) 



so that (5.11) becomes 



u'{ + Ul 



(5.13) 

ay + 13 

Clearly we could regain the equation of a simple harmonic oscillator by 
means of the further change of variable v = ui — uips, where Uipg IS a 
particular solution of (5.13). Consequently the Kepler problem with drag 
differs from the standard Kepler problem in that the second equation of the 
reduced system, viz 

4 = -a, (5.14) 

has a nonzero right side. This could be eliminated by now making the second 
choice mentioned above. Consequently we make a final change of variables 



Vl 
V2 



Ul - 

U2 + a 



Hsm{y — s)ds 
as — P 



(5.15) 



to obtain the same reduced system, ( |4.9| ), as we had for the Kepler problem 
when the introduced the sensible coordinates. Naturally we obtain the same 
set of symmetries as given in (4.10). 

The Lie point symmetries of the reduced system translate into the sym- 
metries 



r2 

r4± 



dt 



dt 



{ar + I')dt 



2t + 



r{ae - (3) 
Idt 



dt + 



r{ae - (3) 
:±^^dt 



ar + /' 

ae-13 

dt+{r + 



r{ae - (3) 



dt + 



ae-(3 



a9-P 

dr 



dr 
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+e 



±2ie 



ar + 1 
r{ae-l3) 
ar + /' 
ae-13 



T 2 + 



r{ae - 13) 



ae-(3 



dr 



do 



{r{ae -(3) + I)[ ""^.^ ^1 T M 3 + 



21 



+ar + /' + ^{aO - B) 

e 
I 



r{ae-p) 
dt\dt + e^'^ {r{ae - f3) + I) 



r{ae-(3), 

ar + /' 



=FM ^ + 



ae-13 



dr-e^'^ [r{ae- l3)+I]dB, 



_al3-p 
(5.16) 



where / stands for the integral introduced in ( ^.ISl a) and /' its derivative 
with respect to 9, for the original system, ( |5.1D. 

in~(|j 



In addition to the symmetries listed in ( ^.16 ) equation (|5.1| ) has the 
point symmetry dt which was the symmetry used for the reduction of order. 
Consequently we can conclude that algebraically the Kepler problem and 
the Kepler problem with drag are identical. 

In the above derivation we have followed a line of development in which 
observation and experience play major roles in reducing the system (|5.4| ) 
to the simplest possible form. The need for both are considerably obviated 
when the interactive Lie symmetry solver devised by Nucci |2j] is used. 
The equations to be solved suggest the appropriate variables since they are 
the characteristics of the partial differential equations to be solved. We 
illustrate this in the case of the variable related to angular momentum with 
the following tableau 



W3 



W4 = Ui,Wi = U2 



Ui 



W2 = y 
dy 



a 



U2 

,2 , o„-,2\„,2„ 



U2 



{U2 + 2u^)u{u2 - IJ, 



ti'4 = — 2 



W5 = Ui,Wi = U2 



Ui 



-a 
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-fiui + uiui + 2uiu2 

U2 = 2 

W5 = wq- ay 

= Ul,Wl = U2 

ui = 

.. ^ {ay-uif{ul + 2ul) - fiu^ 
{ay - uiYu2 

(We have used , and wq to successively define a new variable which is a 
candidate for selection as ui. We are not introducing additional variables.) 

6 The generalisation of the Kepler problem with 
drag 



(6.1) 



Equation (1.4) has, in two dimensions, the two components of the equation 
of motion 

n2 , 1 f 9' -2 
r = rcy+^ 1 — r— ^gr 

\9 rj 
2r\g r, 

We introduce the variables Wi^i = 1,4 as above. Now the system of first 
order equations in these variables is 

Wl = Ws 
W2 = W4, 

2 1/5' 

^3 = WiW4^ + 2 + ~ j "^3 - f^9Wl 

W3W4 / g' 1 \ 
W4 = . (6.2) 

2wi \g wiJ 

As the new independent variable we take again y = W2- The system ( |6.2| ) 
becomes 

dwi W3 , 

— — = W3 = W^Wi 

dy W4 

dws 1 /g' 3\ ,2 "'I 

— — = W1W4 +2 h - W4W 1 - ng — 

dy \g rJ W4 

- ^'-w,-l). (6.3) 



dy 2wi \g 
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When we substitute for W3, the third of ( |6.3| ) is easily integrated to give 



A 



Wl 

9 



(6.4) 



where A is an arbitrary constant of integration. The right-hand side is the 
same function as the characteristic of the paraboUc partial differential equa- 
tion produced when the system ( |6.3| ) is analysed using the code developed 
by Nucci |2j] and this is an appropriate choice for one of the variables. We 
take ui = Wl = r and U2 = 0{r/gY^'^ so that with the elimination of ^3 
from (|6.3| ) we have, after a certain amount of simplification, the system of 
two equations 



U2U1 



I 

U2 



U1U2 + 2 



0. 



11 

U2U'f 
Ul 



U2 



(6.5) 



We observe that ( |6.5| ) is precisely the system ( |3.13| ) and so we immediately 
introduce the new variable v = ^ — u^/ui to obtain the simpler system ( [4.9| ) 
which has the symmetries listed in ( [1.10 ). In terms of the original variables 
these symmetries are 

Xi = 
X2 = 

X3 = 

X4± = 



Xq± 




(6.6) 
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7 An example with an angle-dependent force 



Sen ||2^ obtained conserved quantities similar to those of the Kepler problem 
for the Hamiltonian 



H 



2 , Pe 



r 



a sm 



1(9 -P) 



a/2 



(7.1) 



in which the potential depends upon the azimuthal angle and fi, a and f3 are 
constants. Subsequently Gorringe and Leach Q showed that the equation 
of motion 

r + gr + he = 0, 



where 



u"{e) + ui9) ^ ^v'je) 

^2 ^3/2 



and h 



VjO) 



or 



re + 2re + h = 



(7.2) 
(7.3) 

(7.4) 
(7.5) 



in plane polar coordinates, could be solved for the orbit equation in a manner 
similar to that of the Kepler problem since it also possessed a Laplace- Runge- 
Lenz vector. The only restrictions on the functions U{e) and V{e) are that 
they be differentiable. 

We make the same reduction as in the previous cases to arrive at the 
two-dimensional system 



// 2 /22 22 / 

wiWiW^^ — 2wi w^^ — WiW^^ = Wih — gwi 
wiW4w'4^ + 2w'iw1 = —h, 



(7.6) 
(7.7) 



where wi = r and Wi = 9 as before. 



The Laplace- Runge-Lenz vector for equation (7^) is 



i- X L — C/r 



U' + 2ri/V 



(7. 



where L := r x r is the angular momentum. If we take the two cartesian 
components of J, viz Jx and Jy, and combine them we obtain 

= (r^e'^ -U^±i (-r'^i-e - U' - 2r^/V 
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(wfwl 



'I? 



[/) ± i (^-wjwiws -U' - 2wl^^V 



U\±i 



U 



Wl 



(7.9) 



We see that, when we write the components of the Laplace-Runge-Lenz vec- 
tor in this form, we have the same structure as for the standard Kepler 
problem. (One could call the components the Ermanno-BernouUi constants 
in honour of the original discoverers of these conserved quantities.) Imme- 
diately we have the clue to the identification of one of the new variables and 
we let 

Ul = wfwl -U = U (7.10) 



Wl 



U 

so that the Ermanno-Bernoulli constants for ( |7.^ ) are 

J± = (ni ± iu[) e'^'y. 



The identification of the second variable is more delicate. Equation (|7.7| ) 
can be written in terms of the magnitude of the angular momentum, L, as 



(7.11) 



LL' = -wl^'viy) 
and, when ( [7.10|) is taken into account, this becomes 

Viy) 



L' 

0= ^ + 



L2 {ui + U{y)fl^' 



(7.12) 



(7.13) 



; From ( [tIiD we have 



Ul 



J cosy 



ly 



+ J_e 



-ly 



(7.14) 



since J_ = J\_ and we have written J = [ J+ 1 = | J- [ . Then we can use ( |7.13|) 
to define a new variable 



1 



{J cosy + U{y))^l'^' 
The reduced system of equations is 



(7.15) 



Ul + Ul 



u'2 








(7.16) 
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which is just the reduced system we obtained for the standard Kepler prob- 
lem and so it has the symmetries given in ( 4.10I ). These translate to 



r4± 



r6± 



r^9dtj dt + 2r^9dr 

2U' 3Vr^e 
+ 



r.202 U + Jcosl 



dt 





r Udf 


dt + 


\ u 1 


t - 




r — 




1 r^02 







dt + 

dr 



U' 2Vr^e 

+ 



"26)2 U + Jcosl 



dr - da 



2 I / ^dt dt + 

j,3Q2 I j,2Q2 

U'TiU 



,±2ie 



U + JcosO 
2Vr^^ 



[/ + J cos ( 



± ir + 



2U' 1 



U 

^Iq2 



dtjdt 

dr + e^^'^de 



dt\dt 



2Vr'^e 
U + 7008 6* 



{u'±i(rH^-U)) (l 



u 

r3Q2 



r^d^ - U 



do 



(7.17) 



for the original system, ( |7.2D . In addition there is the symmetry, dt, which 
was used for the reduction of order. 

Again we see the very close connection between the structure of the 
Ermanno-Bernoulli constants and the appropriate variables for the reduction 
of order. 



8 Conclusions and observations 

For the Kepler problem in three dimensions we obtain the same symmetries 



as in (3.14) with the addition of 



X 



10 



du3 ) 



^.11 



where is the azimuthal angle, (j). Consequently our analysis in the lower 
dimensional configurational space is justified by the result that the addi- 
tional dimension simply adds another ignorable coordinate to the original 
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system and so a trivial first order ordinary differential equation to the re- 
duced system. 

In this paper we have examined the process of reduction of order intro- 
duced by Nucci [^] to derive the additional nonlocal symmetries required 
for the complete specification of the Kepler problem in the context not only 
of the Kepler problem but also in some generalisations which have appeared 
in the literature and which possess certain characteristics in common with 
the Kepler problem. In particular we have found that the possession of a 
conserved vector similar to that of the Laplace-Runge-Lenz vector, whether 
or not the magnitude of the angular momentum is conserved, leads in all 
cases to a reduced system consisting of the simple harmonic oscillator and a 
trivial first order ordinary differential equation. In the reduced system the 
Lie point symmetries can be written in a fairly simple fashion. When trans- 
lated to the original system, they are not so simple in appearance. However, 
one can determine the algebraic properties of the several systems studied 
from those of the reduced system provided one adds the Lie point symmetry 
used in the reduction of order, viz dt- The Lie algebra of Lie point symme- 
tries of the reduced system is Ai © sl{3,R) and consequently the original 
systems each have the algebra 2Ai © sl{3,R) since the symmetry used in 
the reduction of order has a zero Lie bracket with the other symmetries. 

Considering the results obtained in this paper we can envisage a reversal 
of the procedure. Instead of taking a system which has a vector of the type 
of a Laplace-Runge-Lenz vector we could simply commence with the reduced 
system and introduce some transformation of the two "reduced" variables 
and an Ansatz on the relationship defining the variable we have been denot- 
ing by W3. One could expect to obtain many "lame ducks"! However, there 
is one aspect which has the potential for some application. Many of the sys- 
tems for which Laplace-Runge-Lenz vectors have been obtained do not have 
a known Hamiltonian representation. By the procedures of transformations 
treated in this paper one could seek to commence with the Hamiltonian of 
the Kepler problem and find Hamiltonians for the other systems. In the case 
of the system ( |7.3| ) the existence of a Hamiltonian has been shown only in 
the restricted case treated by Sen. The attractions for the applications in 
quantum mechanics are obvious. 

In the reduced system the components of the Laplace-Runge-Lenz vector, 
the Ermanno-Bernoulli constants, are simply the two linearly independent 
first order invariants of the simple harmonic oscillator. In fact in the reduced 
system we have a separation in the new variables of the Ermanno-Bernoulli 
constants in the second order equation and the conservation of a generalised 
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angular momentum in the first order equation. We recall that for higher 
dimensional oscillators there exist the conserved components of the Jauch- 
Hill-Pradkin tensor |jl3|, |^] which play an important role in the description 
of the orbit and their time-dependent counterparts which give the actual 
trajectory of the particle [17, 18, 21|. Naturally these tensors have no role 
to play in the type of problem considered in this paper. However, it is 
intriguing to ponder the identity that the corresponding problem of Kepler 
type would have. (For a recent contribution to this more general problem 
see 111.) 
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